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We establish a correlation for the symmetry energy at saturation density So, slope parameter L 
and curvature parameter K synl based on widely different mean field interactions. With the help of 
this correlation and available empirical and theoretical information, the density dependent behavior 
around the saturation density is determined. We compare the results obtained with the present 
approach with those by other analyses. With this obtained density dependent behavior of the 
symmetry energy, the neutron skin thickness of 208 Pb and some properties of neutron stars are 
investigated. In addition, it is found that the expression S(p) = So(p/po) J or S(p) = 12.5 (p/po) + 
C p {p/po) 1 does not reproduce the density dependence of the symmetry energy as predicted by the 
mean-field approach around nuclear saturation density. 

PACS numbers: 21.65.Ef, 21.65. Cd, 26.60.Gj 
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I. INTRODUCTION 

Our knowledge about the features of nuclear matter at 
saturation density po is based primarily on the masses of 
nuclei, like the density po = 0.16 fm -3 , energy per par- 
ticle a v = — 16 MeV and symmetry energy S(p — po) — 
28 — 34 MeV [l|. However, the variation of the symme- 
try energy with density is still intensely debated 
The symmetry energy which characterizes the isospin- 
dependent part of the equation of state (EOS) of asym- 
metric nuclear matter, plays a crucial role in many issues 
of nuclear physics as well as astrophysics. It relates the 
heavy ion reactions [6l-fl3j , stability of superheavy nuclei 
[l4j . fusion cross sections [l5[ and structures, composi- 
tion and cooling of neutron stars fl6T - [l9l ]. Many theo- 
retical and experimental efforts have been performed to 
constrain the density-dependent symmetry energy 0-Q ■ 
The energy per particle of nuclear matter with density 
P = Pa + Pp and asymmetry 6 — (p n ~ p P )/p is usually 
written as e(p,S) = e(p, 0) -I- S(p)5 2 , where p n , p p and p 
are the neutron, proton and nucleon densities. Around 
the nuclear matter saturation density po, the symmetry 
energy S(p) can be expanded to second order in density 
as 
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are the slope and curvature parameters at po that gov- 
ern the density dependence of S(p) around po- So de- 
scribes the symmetry energy at density po Recently, 
some progress has been made in determining the density 
dependence ofS(p) around the saturation density po from 
nuclear isospin diffusion, double n/p ratio in intermedi- 
ate energy heavy-ion collisions, pygmy dipole resonance, 
neutron skin thickness and the nuclear binding energy. 
A description of isospin diffusion data with a symme- 
try energy of S(p) — So(p/po) 7 with 7 = 0.69 — 1.05 has 
been obtained by using an isospin-dependent Boltzmann- 
Uehling-Uhlenbeck (IBUU) transport model [20], and a 
value of 7 = 0.5 is inferred from the preequilibrium 
neutron and proton transverse emissions comparing to 
IBUU transport calculations 21] . The giant dipole res- 
onance (GDR) of 208 Pb analyzed with the Skyrme in- 
teractions implies 7 = 0.5 — 0.65 (22J. The S(p) ex- 
tracted from more than 2000 measured nuclear masses 
gives 7 = 0.6 — 0.8 [23J. More recently, the analysis of 
isospin diffusion and double ratio data involving neutron 
and proton spectra by an improved quantum molecular 
dynamics transport model suggests 7 = 0.4 — 1.05 with 

S(p) = 12.5(p/p ) 2/3 + C p (p/p ) 7 M- Although sig- 
nificant progress was made in determining the symmetry 
energy, it remains an open question nowadays. 

It has been established that the neutron skin thickness 
AR np , given by the difference of neutron and proton root- 
mean-square radii of heavy nuclei, correlates linearly with 
the slope L around the saturation [25 w 30] . Although the 
theoretical predictions on S(p) with the current nuclear 
mean field methods and neutron skin thickness are ex- 
tremely diverse, this correlation is universal in the realm 
of mean field theory as it is based on widely different nu- 
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clear functionals [29|, |30j ■ Based on the similar idea, in 
this work, we try to extract a relation of the three quan- 
tities So, L and K sym in widely different mean field in- 
teractions to constrain the density dependent symmetry 
energy S(p). This work is organized as follows. In Sec. 

II, we extract a relation between So, L and -Ksym univer- 
sally within the mean field framework since it is based on 
widely different nuclear mean-field interactions. In sec. 

III, by employing this relation and other considerations, 
we determine the density dependence of the symmetry 
energy around the saturation density. With the obtained 
density dependent behavior of the symmetry energy, the 
neutron skin thickness of 208 Pb and some properties of 
neutron stars are investigated. Finally, a short summary 
is given in Sec. IV. 
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II. ESTABLISHMENT OF THE RELATION FOR 
THREE QUANTITIES So, L AND K Bym 

Let us first establish a relation of the three quantities. 
With the relation S(p) = So(p/po) 7 descri bing the den- 
sity dependence of the symmetry energy [201 ] . one can 
obtain L = 3So7 and K sym = 9607(7 — 1)j and thus a 
correlation of So, L and K sym can be derived 



FIG. 1: (Color online) A S = S - L/(3 + K sym /L) (rectan- 



gle symbols), A So = So — 12.5 - 
(circle symbols) and ASo = So 



(L — 25) / (3L + Ks- 



50) 



2.6 - L/3 - -?f sym /15 (trian- 



gle symbols) within the density dependent behavior S(p) = 

S (p/po) 7 , S(p) = 12.5 (p/po) 2/3 + C p ( P / Po y and DDM3Y 
shape, respectively. The horizontal ordinate denotes the se- 
quence number for the interactions mentioned in text. 
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For the other density dependent behavior S(p) — 
12.5 (p/p ) 2/3 + C p {p/poV, we have L = 25 + 3(S - 



12.5)7 and K s 



-25 + 9C P 7(7 - 1). Then, the corre- 



lation of So, L and K sym takes the form 



So = 12.5 
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The shape of the density dependence of the symme- 
try energy S(p) that from the density- dependent M3Y 
(DDM3Y) interaction 32] can be written as S(p) = 

C k (p/p ) 2/3 + C 1 (p/ Po ) + C 2 (p/pof 3 with C fc = (2 2 / 3 - 

i)!^ = 13 -° MeV - 5 o, L and K sym can be expressed 
by parameters C\ and Ci with So= 13.0 + C1+C2, L = 
26.0 + 3Ci + 5C 2 and K sym = -26 + 10C 2 . Therefore, 
the correlation is given by 
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We now test whether Eqs. (0), ([3]) and flU work well 
or not by using widely different mean field functionals in- 
cluding relativistic and non-relativistic versions. As done 
in Ref. [3ll j . to prevent eventual biases in our study, 
we avoid including more than two models of the same 
kind fitted by the same group. We also avoid models 
yielding a charge radius of 208 Pb away from experiment 
data by more than 1% as in Ref. 31] since we will 
study the neutron skin thickness of 208 Pb. The inter- 
actions we used here are (1)LNS1, (2)LNS5, (3)MSL0, 



(4)SIV, (5)SkT4, (6)T6, (7)SkP, (8)SkM*, (9)SkX, 
(10)PK1, (11)D1S, (12)SLy4, (13)FSUGold, (14)SkMP, 
(15)SkI5, (16)NLSH, (17)TM1, (18)NL3, (19)NL1, 
(20)Sk255, (21)DDME1, (22)DDME2, (23)DDM3Y, 
(24)PC-F1, (25)Ska, (26)SV, (27)QMC, (28)MSkA, 
(29)SkI2, (30)MSk7, (31)HFB-17, (32)BSk8, (33)BSkl7, 
(34)GM1, (35)GM3, (36)Sk272, (37)v090. The root- 
mean-square deviations given by Eqs. ([2|), ((3|) and (j4]) 
are 38.65 MeV, 65.29 MeV and 2.94 MeV for S value, 
respectively. In order to give a clearer show, we define 
a quantity ASo to describe the differences between the 
left hand side and the right hand side of Eqs. ©, §5§ 
and and plot them in Fig. 1. As can be seen, 
ASo with Eq. ^ approximates zero for these widely 
different interactions while ASo with Eqs. ((2|), ([3]) tends 
to deviate from zero considerably. Therefore, Eqs. 
and ^ can not be taken as accurate expressions for 
the description of the correlation of So, L and K sym 
displayed by the mean-field calculations. It is noted 
that the minimum value of K sym is — 9So/4 m — 72MeV 
with S(p) = S (p/po) 7 and -25 - 9C p /4 « -69 MeV 

with S(p) = 12.5 {p/p ) 2/3 + C p (p/po) 1 - Many inter- 
actions, however, provide K synl < —100 MeV (as will 
be seen in Fig. 2), which leads to large discrepancies 
of Eqs. @ and ([3]) when compared to the mean-field 
predictions. This fact also indicates S(p) — So(p//9o) 7 
or S(p) — 12.5 (p/po) 2 ^ 3 + C p (p/po) 1 is not suitable to 
describe the density dependent behavior of the symme- 
try energy around po as predicted by the mean field ap- 
proach. In particular, they cannot describe a very soft 
symmetry energy due to their monotonous increase with 
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the density. However, Eq. (HJ is much better to be taken 
as a relationship for the correlation of So, L and -K sym . 
Noting that the S(p) with the DDM3Y shape is still ap- 
proximate, Eq. Q can be further improved to obtain the 
least deviation. Replacing the index 5/3 in the DDM3Y 
shape by a coefficient 7, one has So = a + L/3 + bK sym - 
By performing a least-squares fit with the calculated So, 
L and K sym using the interactions above, the values of 
the parameters are a — 3.9199 and b = —0.07323 with a 
rms deviation of 2.12 MeV for So value which is slightly 
different from that within the DDM3Y shape, and hence 
the relation is give by 

L = -11.76 + 35 + ^ff, (5) 
4.55 

with a rms deviation of 6.35 MeV for L value. This for- 
mula can be considered as a universal one within the 
mean-field framework since it is based on widely differ- 
ent nuclear mean-field interactions. As a consequence, 
the analytical and simple expression reported in Eq. ([5]) 
describes in very good approximation the high correla- 
tions displayed between So, L and K sym arising from the 
predictions of the representative set of employed mean- 
field models. The corresponding expression for the den- 
sity dependence of the symmetry energy is given by 

S{ P ) = C k (^) 2 '\c x U) + cJ^] m , (6) 
\PoJ \PaJ \PoJ 

where Ck = 17.47 is larger than that in the DDM3Y 
shape, which perhaps can be regarded as the mass in the 
kinetic energy being replaced by the effective mass m*. 
In this step, parameters C\ and C2 remain unknown. 

III. DENSITY DEPENDENCE OF THE 
SYMMETRY ENERGY AROUND p WITHIN 
THE OBTAINED RELATIONS 



TABLE I: The So values obtained from various independent 
studies in other references. 



Reference 


So (MeV) 


Reference 


So (MeV) 


Ref. [3, 20J 


31.6 


Ref. [23J 


31.1 ± 1.7 


Ref. [24] 


30.1 


Ref. [33] 


32.4 ± 1.1 


Ref. [34] 


32.0 ± 1.8 


Ref. [35] 


32.0 


Ref. [36] 


30.048 ± 0.004 


Ref. [4J 


32.3 ± 1.3 



Of the three quantities So, L and -Ksym, So value is 
relatively well-known. Table I lists the recent So value 
coming from recent various studies. The largest range 
of So = 31.6 ± 2.2 MeV from Table I will be used in 
the present study. Yet, we have to call for an additional 
condition to constrain the detailed L and -Ksym values. 
Recently, Centelles et at found that the symmetry en- 
ergy (coefficient) a sym (^l) of a finite nucleus with mass 
number A is approximately equal to the symmetry en- 
ergy S(pa) of nuclear matter at a reference density pa, 
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FIG. 2: Correlation of the curvature K sy m with the slope L of 
the symmetry energy at po- The line gives the fitting result 
with the correlation coefficient r — 0.972. The inset displays 
the correlation between the neutron skin thickness Ar np in 
208 Pb and K Byul value. The line gives the fitting results with 
r = 0.945. 



namely, S(pa) = a sym (^4) [2!|- For a given nucleus, the 
Pa is determined, such as pa = 0.1 fm~ 3 for 208 Pb. Here, 
this relationship is employed to investigate the correla- 
tion between K sym and L. a sym (^4) of a finite nucleus is 
given by 

aSym(A)= l + ^-V3 ' WithK -lt' 

where Q is the surface stiffness that measures the resis- 
tance of the nucleus against separation of neutrons from 
protons to form a neutron skin. In Ref. [30l |. it is shown 
that So/Q ratio displays a linear relationship with L. 
Then, parameter k can be written as k = mL + n. By 
combining Eq. ([7]) and the expression of the symmetry 
energy obtained with the DDM3Y shape aforesaid, the 
relationship of S(pa) = a sym (A) can be converted into 
L 2 + aiL-Kgym + a2-K sym + a^L + 04 = with new coef- 
ficients ai,d2,a3 and 04 by replacing So, C\ and C2 by 
L and -K sym . We fit the results from the nuclear mean 
field calculations with the interactions mentioned above. 
However, in the fitting process, it is found that the first 
two terms can be neglected in the region under consider- 
ation. No matter whether or not the first two terms are 
taken into account, one can obtain the same root- mean- 
square deviation v< a 2 > = 29.8 MeV for K sym value. 
Accordingly, this correlation can be further simplified as 
a linear relation, as visibly shown in Fig. 2. By per- 
forming a two parameter fitting to this relationship, we 
have 

K sym = -307.862 + 3.292L, (8) 



with V < c 2 > = 29.8 MeV for K sym value and the cor- 
relation coefficient r = 0.972. The results of fitting are 
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presented in Fig. 2 by the line. One readily sees that the 
predicted K sym with the mean field approaches varies 
largely from -300 MeV to 150 MeV. By combining Eq. 
© and Eq. © together with S , = 31.6±2.2 MeV, in 
the present work, some information on the density de- 
pendence of S(p) is investigated. 

Fig. 3 illustrates the present estimated L values com- 
pared with those from other independent approaches. 
One can find that the present finding has a remark- 
able overlap with but on the whole slightly softer than 
the very recent results of Carbone et al. 01, Liu et 
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al. [23j and Tsang et al. 24]. This, to a large ex- 
tent, perhaps stems from the relation of Eq. ([5]) based 
on the formula of S(p) that can characterize the very 
soft symmetry energy in a much more reasonable man- 
ner. In the other analysis with S(p) — Sq(p/ po)* 1 or 
S(p) = 12.5 (p/po) 2 ^ 3 + C p (p/po) 1 , the slope parameter 
of the symmetry energy at saturation density should be 
overestimated. To give a clearer explanation, we present 
an example here to show the overestimation of the L val- 
ues with S(p) = So(p/ pq) 1 . In Ref. (33|, Danielewicz 
and Lee extracted the mass dependent symmetry energy 
coefficients of finite nuclei a sym (A) = Sn (l + k^4 -1 / 3 ) 
with S = 27.39 MeV and k = 1.28 (37]]. With the rela- 
tion a sym (A) = S(pa) proposed in Ref. [29j |. if the den- 
sity dependent behavior S(p) = So{p/ po) 1 is adopted, 
one obtains L = 34 MeV. However, if the DDM3Y shape 
or the modified one (Eq. ©) is applied, we obtain L = 21 
MeV and L — 24 MeV, respectively, lower than that 
stemming from S(p) = S {p/poV- With S = 29.4-33.8 
MeV in Table I, the L and AT sym values obtained with the 
present method are 56 ± 24 MeV and -125 ± 79 MeV, 
respectively. It is interesting to see that the window re- 
ported in reference of Warda et al. [30] for the analy- 
sis of anti-protonic atoms data and that in the present 
manuscript are exactly the same, as shown in Fig. 3. 
The So value with less uncertainty will lead to narrower 
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FIG. 3: Comparison between the L values obtained in the 
present work and those from other recently various analyses 
including Carbone et al. [H, Liu et al. l23ll. Tsang et al. [24|l. 
Warda et al. \3uh, D anielewicz and Lee [33], Shetty et al. [3| 
and Chen et al. 12011. 



spectively. The result of the fitting is presented in the 
inset of Fig. 2 by the line. A larger K sym value im- 
plies a thicker neutron skin. The linear relation allows 



one to extract K s . 



with the measured AR np . There- 



fore, once an accurate measurement of AR np is achieved, 
L as well as K sym values can be determined simultane- 
ously. As a consequence, a richer information about the 
density dependent behavior of the symmetry energy can 



be achieved. Here K Byir 
the AR nr> . With the K s 



is used in turn to determine 
m values obtained above with 

208t 



S = 31.6 ± 2.2 MeV, the AR np in 2m Pb is estimated 
windows of L and K sym values within our approach. The to be 0.185 ± 0.035 fm^With the relationship between 
errors for L and K sym arising from the errors of the pa- 
rameters in the fitting are less important compared with 
these resulting from the uncertainty of So value. There- 
fore, for simplicity and clarity, we only selected the op- 
timal fitting results. The small deviation and large cor- 
relation coefficient indicate that the parameters are well 
constrained by the information used in the fitting proce- 
dure. 

It is well known that L shows a linear relationship with 
neutron skin thickness AR np of finite nuclei. In Ref. 
[3lj . this linear correlation is given by AR np = 0.101 + 
0.00147L for 208 Pb. Because of Eq. ©, K sym should also 
show a linear correlation with the neutron skin thickness 



AR 



Using the aforesaid interactions for finite nuclei, 
we obtain the correlation between AR np in 208 Pb and 
-Ksym by performing a fitting procedure 

AR np = (0.236210 ±0.003991) 

+ (0.000415 ± 0.000026)if sym , (9) 

where AR np and K sym are in units of fm and MeV, re- 



AR np and L in Ref. [3l| together with the L values of 
56 ±24 MeV that we obtained above, the AR np for 208 Pb 
is 0.183 ± 0.035 fm, which is quite consistent with that 
from AR np — K sym correlation. This fact indicates the 
justification of our approach to a large extent. Again, 
our calculated AR np for 208 Pb is in good agreement with 
the value from analysis of PDR 0, S] • The neutron-rich 
skin of a heavy nucleus is related to the properties of 
neutron star crusts. For instance, the thicker the neu- 
tron skin is, the thinner the solid crust of a neutron star 
[l6l ] . Horowitz and Piekarewicz proposed that if AR np for 
208 Pb is greater than about 0.24 fm, the electron fraction 
becomes large enough to allow the direct URCA process 
to cool down a 1.4M Q neutron star [l6|]. Our calculated 
ARnp is too small to allow the direct URCA process in 
this canonical neutron star to occur. An almost linear re- 
lationship between the AR np and the critical density p c 
of a phase transition from nonuniform to uniform neu- 
tron rich matter is put forward in Ref. flol ] with an 
approximate relation of p c w 0.16 — 0.39Ai?„ p ( 208 Pb). 
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Using this relation combining with AR np in 208 Pb de- 
duced above, one obtains p c s» 0.09 fm~ 3 which is in 
accord with p c w 0.096 fm -3 of the microscopic EOS of 
Friedman and Pandharipande [38| . These two examples 
display the applications of the correlations we employed 
in understanding the physics of compact objects. 
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uncertainty. As suggested by Daoutidis and Goriely [4(j , 
the PDR strength measurements nowadays cannot yield 
quantitative insight on the symmetry energy, and thus 
both theoretical approaches as well as phenomenological 
interactions need to be further improved to investigate 
the symmetry energy. Some authors even proposed that 
the PDR strength is weakly correlated with the neutron 
skin thickness of heavy nuclei [4l[ and hence the density 
dependence of the symmetry energy. Finally, we would 
note that compared with the transport models and the 
PDR measurements, our approach is much more straight- 
forward. Yet, our approach is based only on mean field 
models and could only gain some limited information 
about the symmetry energy. Opposite to our method, 
the analysis of the experimental data within microscopic 
techniques, such as the random phase approximation for 
the study of the PDR, may allow one to investigate the 
internal structure and dynamics of the nucleus, and hence 
one can obtain more detailed knowledge. These different 
approaches can validate or complement each other to get 
more compelling results. 



FIG. 4: (Color online) Behavior of the density dependent 
symmetry energy with different So values from Eqs. {lj (dash 
lines) and (0 (solid lines), respectively. 

When So value is fixed, the parameters C\ and Ci 
in Eq. ^ are accordingly determined, and hence the 
density dependence of the symmetry energy in the form 
of a modified DDM3Y shape is obtained. Within this 
density dependent behavior of S(p), one may investigate 
the symmetry energy at subsaturation densities. Fig. 4 
presents the symmetry energy S(p) as a function of the 
density with several detailed Sq values. The solid (dash) 
curves denote the calculations with Eq. ([6]) (Eq. (JTJ). As 
can be seen, S(p) shows different behaviors with different 
Sq values. The symmetry energy tends to become stiffcr 
with a larger 5*0 value. When 5*0 is selected to be 33.8 
MeV, the relation of S(p) — p is almost linear. This ex- 
plains the fact that the relativistic mean field models have 
traditionally predicted a stiff symmetry energy-these in- 
teractions give larger So values. The density dependent 
behavior stemming from Eq. ((5]) and Taylor expansion of 
Eq. (jTJ) mostly coincide with each other, which indicates 
that the Taylor expansion of Eq. ([1} can be applied in a 
wide density range when Eq. (|6|) is employed for the de- 
scription of S(p). In addition, consistent with the mean- 
field results of Refs. [Ill [27|, [2!| , ^ is found that different 
curves give nearly the same value for the symmetry en- 
ergy at a density around p = 0.1 fm -3 . The obtained 
symmetry energy S(p) = 23.6 ± 0.02 MeV at a nucleon 
density of p = 0.1 fm" 3 agrees with 23.3 - 24.9 MeV by 
an analysis of the GDR of 208 Pb with the Skyrme inter- 
actions [22[ but is larger than the values of 21.2 — 22.5 
MeV obtained with analysis of the GDR of 132 Sn within 
relativistic mean field models [39| . Moreover, compared 
with that from the PDR, our result includes a much less 



IV. SUMMARY 

Based on the similar idea in Ref. [25M30l ] that the neu- 
tron skin thickness correlates with the slope parameter 
L around the saturation density in mean field models, 
we have established a relation for three quantities So, 
L and K sym in widely different mean field interactions. 
With this relation and other constraint conditions, the 
density dependence of the nuclear symmetry energy S(p) 
has been investigated in the present work. With the ob- 
tained density dependence of the symmetry energy, the 
neutron skin thickness of 208 Pb and some properties of 
neutron stars were analyzed. The main conclusions are 
as follows, (i) It is not suitable to take the form of 

S(p) = So(p/poV or S(p) = 12.5 (p/po f + C p (p/po) 1 
to describe the behavior of the symmetry energy as pre- 
dicted by the mean-field approach around the nuclear 
saturation density while the shape from the DDM3Y is 
much better. (ii)Based on the latter formula, we have 
provided an analytical and simple expression for the high 
correlation shown by 5*0, L and K sym within the mean 
field framework. Therefore, if the symmetry energy at 
saturation density po is known, the values of AT S ym and L 
can be unambiguously related by using Eq. ([5]). The L 
and K sym values in the present study are 56±24 MeV and 
-125 ± 79 MeV respectively with S = 31.6 ± 2.2 MeV 
yielded in other references. (iii) The neutron skin thick- 
ness of 208 Pb displays a linear correlation with K sym . 
Thus, once the neutron skin thickness is measured ac- 
curately, not only slope parameter L but also curvature 
parameter K sym of the symmetry energy around the sat- 
uration can be determined. Thus, a richer information 
about the density dependent behavior of the symmetry 
energy can be achieved. AR np of 208 Pb is 0.185 ± 0.035 
fm with the -Ksym value in item (ii) which is too small 
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to allow the direct URCA process in 1.4M neutron star 
to occur. (iv) Within our approach, the behavior of S(p) 
against p around the normal density almost relies on the 
symmetry energy at saturation density po- We show 
that a large So value leads to a stiff symmetry energy. 
They, however, naturally provide almost the same result 
of about 23.6 MeV at p = 0.1 fm" 3 . 
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